Nondeterminstic ultrafast ground state cooling of a mechanical resonator 
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We present an ultrafast feasible scheme for ground state cooling of a mechanical resonator via repeated 
random time-interval measurements on an auxiliary flux qubit. We find that the ground state cooling can be 
achieved with several such measurements. The cooling efficiency hardly depends on the time-intervals between 
any two consecutive measurements. The scheme is also robust against environmental noises. 
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Introduction.- In the quantum regime, ground state coohng 
of a small thermal object is an intriguing challenge and one 
of the most desirable quantum technologies. Physically, the 
cooling process can be formulated as a transformation from 
a thermal state of the small object into its ground state. The 
transformation is irreversible and cannot be performed when 
the object is isolated. 

A mechanical resonator (MR) is a small mecroscopic me- 
chanical object, usually coupled to an auxiliary setup. The 
physical realization of its quantum ground state has be- 
come more and more important in ultrahigh-precision mea- 
surements, classical to quantum transitions, preparations of 
non-classical states, quantum information processing |[ll-l3|- 
Throughout the years, considerable number of optomechan- 
ical i4l 4lOll and electromechanical ll ll4l6ll models has been 
proposed for achieving their ground state cooling. Examples 
are a bang -bang cooling through a Cooper pair box 0121 . a 
single-shot state-swapping cooling via the superconductor via 
a superconductor ifisll . Recently, some of us ifioll proposed 
a ground state cooling scheme of MR in an optomechanical 
system by controlling fast the optical drives. The best studied 
ground state cooling protocol is the sideband cooling 
IgD designed originally for cooling the atomic spatial motion. 
This cooling is now widely used for the MR cooling exper- 
iments and the recent record for the mean phonon number 

IS 3.8 in 

Theoretically, a MR could be cooled down to its 
ground state in the resolved sideband limit, with the mean 
phonon number less than 1. Since ground state cooling is 
not yet achieved experimentally and seems to become much 
harder when close to ground state, new cooling approaches 
remain desired. 

This work presents a new protocol for MR cooling using re- 
peated projective measurements on an auxiliary qubit. While 
it may be used in any stage of the cooling process of a MR, 
our protocol is aimed at the ground state cooling. Our study 
starts with the cooling with repeated equal time-interval mea- 
surements on the qubit, introduced for purification of quan- 
tum states ITtIi or measurement-based entanglement genera- 
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FIG. 1: (Color online) Schematic diagram of the circuit. The top 
left part (red) is the coupled flux qubit-MR system, where each cross 
denotes a Josephson junction and the bar denotes a doubly-clamped 
MR. The interaction between MR and flux qubit is modulated by an 
in-plane magnetic field Bq. The top right part (blue) is a Josephson 
bifurcation amplifier (JBA) formed by a dc SQUID shunted by a ca- 
pacitance. The bottom part (green) is bias to control the qubit energy 
gap as well as the coupling strength between the JBA and the qubit. 



tion iflil [19I1 . While equal time-interval repeated measure- 
ments work well for ground state cooling, we find, unexpect- 
edly, that the cooling efficiency is even much better when the 
repeated measurements are taken randomly. The ground state 
can be reached in several such measurements in very short 
time. We give explanation on the unexpected phenomenon. 
These results suggest that our protocol is completely robust 
against measurement operational eiTors. In addition to this 
great advantage, the scheme is also robust against environ- 
mental noises. 

Model - We employ a gradiometer-type flux qubit 11261 127 1 



as our auxiliary qubit, though our scheme may be applicable 
to any two-level system coupled to a MR. Fig.[T]is a schematic 
diagram of our cooling setup. The doubly-clamped MR (the 
red bar in Fig. [T), is embedded in a flux qubit circuit which 
composed of three superconducting loops with four Joseph- 
son junctions (JJs). An in-plane magnetic field Bq induces 
qubit-MR coupling via Lorentz force OOll . The top right blue 
(bottom green) part is to measure (operate) the qubit as de- 
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scribed later. 

The Hamiltonian of the flux qubit can be written as Hg = 
HAax/'i + Heaz/2, where and are the Pauli matrices 
in the basis of two persistent current states |t) and |^). Here, 
hA is the tunnehng ampHtude between the two states. The 
bias energy he Hnearly depends on external flux bias and in 
our case is set to zero by pre-trapping one flux quantum $o in 
the loop lUF 
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The MR is modeled as a single-mode harmonic oscillator 
with a high-Q mode of frequency u!„i and effective mass m. 
The entire system is characterized by the Hamiltonian [|2^ 



H = Hq + tajj-ino) a — hgia + a))az^ 



(1) 



where a and are annihilation and creation operators for the 
MR with frequency LOm- The last term denotes the interaction 
between the MR and the flux qubit. The coupling constant is 
g = BoIpLo with Bq the magnitude of the in-plane magnetic 
field, Ip the magnitude of the persistent current in the loop, 
and Lo flie length of the MR. 

Here we consider a MR with fundamental mode frequency 
ujm ~ 27r X 100 MHz. The qubit is tuned into resonance 
or near resonance with the MR by monitoring the tunneling 
A ini 12811 . The coupling constant g, e.g., 27r x 1 MHz, is 
much smaller than the qubit frequency such that the rotating 
wave approximation (RWA) can be used to reduce the Hamil- 
tonian Eq. ([T]i to the standard Jaynes-Cummings (JC) Hamil- 
tonian 

H ~ huj„ia) a H — ^ ^9{0'<^+ + (2) 

where dz,± are the Pauli operators in the new basis of ground 
and excited states: \g/e) = (||) + / - |t))/V2- 

The operator Nc = 0,^0, + \e) (e| in the JC model is con- 
served, such that H can be represented by the direct sum of a 
one-dimensional block in the basis \0,g) and two-dimensional 
submatrices in pairs of bases \n, g) and \n — 1, e) when n > 
1. The Hamiltonian ^ can be therefore diagonalized (we set 

h = i) 



H = -^\0,g) {0,g\ 



EE 

n>l s=± 



£„ ns) {ns 



(3) 



where the dressed eigenstates |n+) = cos6'„ |n — l,e) + 
sm6n\n,g), |n— ) = sin0„|n— l,e) — cos 9n\n,g) and 
the corresponding eigenvalues are = {n — l/2)uJm ± 
y/{A — a;„i)2/4 + g^n. Here 6n satisfies the equation 
tan 26n — 2gy/n/ (A — u!m) for n > 1. 

Unitary evolution and repeated measurements.- We pre- 
pare the whole system initially in a separable state, po = 
\g) {g\ (g) pm- Here pm is the thermal state of the MR (it could 
also be an arbitrary state, where the MR ground state |0) is in- 
cluded). We then perform repeated but unequal time-interval 
measurements (UTIMs) on the flux qubit. The whole system 
evolves under the JC Hamiltonian (|2]i in between the measure- 
ments. The j-th measurement takes place at the time instant 
tj = jr + 5tj, where r is a given time interval and (5<j's 



are random variations in time in the interval (— r/2,T/2). 
The time interval between the (j — l)-th and j-th measure- 



ments is unequal, t. 



5tj — Stj^i (for j > 1 and 



to = StQ = 0). When all Stj = 0, the repeated process is 
reduced into equal time-interval measurements (ETIMs) for- 
mulated in MlTll . After N such ETIMs on the qubits and if 
all measurement outcomes are | g), the density matrix of the 
MR becomes p^^^ (TV) = y/'(r)p„ V;t^(r)/pj^^ (TV), where 
P^''\n) = T[[Vg^{T)p,nV^^iT)] is the survival probabil- 
ity llnllll. Here K,(t) = {g\ er'"''/''' \g) is an effective 
evolution operator only acting on the MR. 

Vg{T) for model ^ is diagonal in the basis {[ti)}: Vg{T) = 
Sn>o ^n{T) \n) {n\. Here the eigenvalues are Ao = e"^'^/2 



and A„ = e 



-l/2)u 



"(cos 17„ 



i sill ri„ T cos 29„ ) with 



= ■\/(A — CL>m)2/4 + g'^n for n > 1. By carefully select- 
ing the time interval r such that cos^ ri„T ^ (for n > 1), all 
the values 



|A„(t)| = \/l — COS^ f2„T COS^ 20n 



(4) 



can be made less than 1, while |Ao(t)| always equals to 1. 

In the large N limit for our specific JC model, the UTIMs 
result in 

Vgin)Vg{r2)...Vg{r^) = I^-I'"^ I") (^1 |0) (0| ' 

where |A„|2 = (|A„(ti)|' |A„(t2)|' ... |A„(t^)|') '^"^ < 1 

for n > 1 and |Ao| = 1. All the density matrix elements of 
the MR will vanish, except that of ground state |0): 

.(^)(7V) = f^nl'^'p'^^ \n) {n\ /P^^\N) (5) 

n>0 



|0) (0| , for ^ 00, 



where = (7t,| pm \n) and the survival probability 

pI"\n) = En^ol^^P'^P™^ ^ Forfliecaseof 
ETIMs, |A„| = |A„(r)|. Repeated measurements drive the 
evolution from the initial state (e.g., the thermal state) to the 
ground state of a MR. The process has the same features as 
those of known ground state cooling methods and is a new 
method for ground state cooling of a MR. Furthermore, dis- 
tinct from usual dynamical cooling schemes, the new cooling 
method is based on repeated non-dynamical measurements. 

Numerical results and robustness.- Consider a 27r x 100 
MHz nano-mechanical resonator with quality factor Qm = 
10^ (7m/27r = 200 Hz), coupled to a flux qubit with the tun- 
neling splitting A ~ LOrn- The MR is initially at its ther- 
mal equilibrium state at the ambient temperature T = 20 
mK, and the corresponding mean phonon number is n(0) = 
l/[exp(fta;,„/fcBT) - 1] = 3.69. Fig. |2] shows for ETIMs 
the mean phonon number n{N), the survival probability 



Pt\N), and the fidelity F!j''\N) = 



\p^^\n)\G) as a 



function of N , the number of measurements. The fines with 
red triangles in Fig. |2]are for the on-resonant A = Wm case 
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decays faster since | A„ | < A„ for 



FIG. 2: (Color online) (a) The average phonon number n{N) after 
equal time-interval measurements (ETIMs) for the initial phonon 
number n(0) = 3.69. (b) The corresponding survival probability 
P^''\N) and fidelity F^''\N) (inset). The gray square lines denote 
the resonant case A = uom- The red triangle lines denote the non- 
resonant case A = l.la;„i. Here g = 0.04ciJm and r = lO/tJm. 



and the lines with gray squares correspond to the off -resonant 
case, A = l.la;„i- The ground state cooHng can be reached in 
both cases. Fig. |2] shows that while the ground state cooling 
requires 60 measurements {n{N ~ 60) w lO^*^), the mean 
phonon number decreases 90 percent with only five measure- 
ments. 
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FIG. 3: (Color online) The average phonon number n{N) (a), the 
survival probability (b), and the fidelity [inset in (b)], after A'^ mea- 
surements in the resonant case A = Um- The lines with red triangles 
denote the equal time-interval measurements (ETIMs) (r = 8/cJm) 
and the lines with blue circles denote the unequal time-interval mea- 
surements (UTIMs) j2^ . Here the initial phonon number n(0) — 
7.84, and the coupling strength g = Q.QAujm- 



The ideal ETIMs, with time interval t, may be difficult 
technically. Fig.[3]shows, for UTIMs, the same physical quan- 
tities as Fig. |2]for the A = cj„i case, but with the higher bath 
temperature T = 40 mK and randomly selected time intervals 
Tj. The mean phonon number is n(0) = 7.84 initially. It is 
noticeable that the ground state cooling can be achieved much 
more efficiently than the ETIMs. It is a remarkable advan- 
tage for experimentalists to achieve the ground state cooling 
of MRs with random time intervals. 

The physical reason is clear In general, for a fixed n, the 

smaller the maximal value A„ = max{|A„(Tj)|}j=i is, 

the faster the term |A„p^ \n) {n\ decays with N . For ETIMs, 
it is unavoidable that there exists A„ {n > 1) very close to one, 
since for a given finite time interval r the periodical function 
cos(ri„r) versus 0,i runs across several times. The cor- 
responding component \n) {n\ therefore decays very slowly. 
Specially, the component |n) (n| will not decay with for 
ETIMs when A„ = 1. However for random UTIMs, the cor- 



responding \Xr, 
any n. 

Our cooling scheme is completely robust against the mea- 
surement operational errors or randomness. It is also robust 
against the relaxation effect of the MR. Our UTIMs can cool 
down a MR, initially at a thermal state with mean phonon 
number < 10 and a quality factor Q„i = 10^, in the time 
interval lOr « 100/wm, with 10 measurements. It is 100 
times less than the MR's relaxation time, ^ l/[n(0)7„i] « 
10000/a;„i. However, since the final survival probability is 
proportional to the initial population probability at the ground 
state, our UTIMs is more suitable for further cooling based on 
a pre-cooled MR at a thermal (-like) state with a smaller mean 
phonon number ^ 10, e.g., by other cooling methods such as 
the sideband cooling. 

We should comment that there is an equal time-interval 
measurement-based cooling of MR proposed in ll29ll using a 
Cooper pair box as the auxiliary, where the measurement ef- 
fect is averaged out and there is no explicit analytical expres- 
sion for the process. However, we find that the method fails 
to achieve the MR cooling for our model after considerable 
number of equal time-interval measurements. We should also 
remark that in the well-known sideband cooling the MR is 
coupled to a high-frequency auxiliary with the faster damping 
rate. The energy flows from the MR to the auxiliary and is 
then lost to the bath quickly. On the contrary, the frequencies 
or the damping rates of the MR and the auxiliary qubit are of 
the same order in our model. Our scheme is non-deterministic 
and based on repeated non-dynamical measurements. 

Implementation of the projective measurements on the flux 
qubit- Our UTIMs on flux qubit can be implemented by 
Josephson bifurcation amplifier (JBA) in a fast and non- 
destructive way. As shown in Fig. [T] a JBA |30:] (blue part), 
is coupled to the flux qubit inductively as the measurement 
device. The JBA consists of a dc SQUID shunted by a capac- 
itance C, subject to a microwave drive Irf cos{ujdt + (Pa)- 
The JBA SQUID loop contains two Josephson junctions of 
identical critical current Iaq and different phase differences 
'PAi, 'PA2 respectively. The current in the loop is I a = 
Ia{^a)cos(pa, with the flux bias in the JBA SQUID, 
set by external coils, Ia{^a) = 2/^o sin(7r$A/3'o)- and 
= {fAi + </'yi2)/2. The JBA circuit forms a driven res- 
onator with nonlinear Josephson inductance. 

Since the JBA is positioned symmetrically with respect 
to the qubit loops 1 and 2, the two loops are coupled to 
the JBA with equal mutual inductance, AIi = il/2. Due 
to the gradiometer design, the total qubit flux is decou- 



pled from the JBA IM |27|]. However, the JBA still cou- 
ples to the qubit loop 3 through its influence on $3. If 
ttM^^a) <C <&0i this influence can be approximated as a 



linear coupling to the dz operator of the flux qubit 0311 13211 . 
that is, an extra interaction term between the qubit and the 
measurement device Hj ~ X{^7ib)(^z cosipA- Here the cou- 
pling coefficient is A(<l>3h) = -(7rM3/A/<I>o)'«('&3b) with 
K($3t) 2aosin(7r($3b/$o)) (dA/da)!^^-, and a = 
2ao cos(7r$3f,/<f>o); "0 is the ratio between the Josephson en- 
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ergy of the smaller junctions and that of the two bigger junc- 
tions in the flux qubit; $3f, is the total flux bias of the loop 
3. Thus an external on-chip bias current 7^3 (green part in 
Fig.[T]) can be used to monitor the coupling strength A($3(,). 

Under a strong microwave drive, the Josephson energy of 
the junction —Eja cos i^a is expanded beyond the harmonic 
approximation and the classical dynamics can be described 
by a Duffing oscillator 13311 . For a certain range of drive 
conditions, the nonlinear oscillator exhibits bistable behavior 
with hysteresis llsol. 331. The two possible stable states cor- 



respond to different oscillation amplitudes and phases, which 
can be distinguished by transmitted or reflected microwave 
signals 134143611 . Switching between the two stable states 
happens when the driving power reaches a threshold. The 
switching probability depends on the value of the nonlinear 
inductance, which in our case depends on the states of the 
qubit. This is because the effective Josephson energy of the 
junctions of the JBA is modified by the interaction Hj as 
EjAio'z) = Ta'^o/Stt — A(<i>3f,)(j2. Therefore, by measur- 
ing the phase of the transmitted microwave signal, the state 
of the qubit is collapsed to one of the eigenstates of its free 
Hamiltonian, in our case. 

In a realistic UTIMs, the measurement time, in the order of 
^ 10 ns, should be considered. This is not a problem when 
measurements are performed as follows: First, the qubit- 
MR interaction is switched off through the in-plane magnetic 
field Bq; then a measurement pulse with readout and latching 
plateau is sent to the JBA to readout the qubit state and induce 
the projection to either \g) or |e) state; after this projection, 
the qubit-MR interaction is switched on again. The process is 
repeated until the MR reaches its ground state. 

Conclusion.- We propose an ultrafast feasible scheme to 
cool a MR to its ground state via repeated random-time- 
interval projective measurements on an auxiliary flux qubit. 
The measurement scheme is almost independent of the initial 
state of the MR. It works when the MR couples with the qubit 
either on-resonance or off -resonance, and even when the cou- 
pling g{t) is time-dependent. The cooling process is robust 
since it can be accomplished in a much shorter time than the 
relaxation time. 

Our scheme significantly simplifies experimental con- 
straints since there is no requirements for the control on the 
time intervals of measurements. In principle, the MR can 
be cooled to arbitrarily low temperature with arbitrarily s mall 
mean phonon number, within very short time. 
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